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SUMMARY 

This  report  discusses  the  characteristics  of  the  gaseous  trail  remain¬ 
ing  behind  a  body  moving  through  the  atmosphere  at  hypersonic  speeds.  The 
purpose  of  this  study  of  the  hypersonic  trail  is  the  development  of  means 
for  ascertaining  those  variables  that  can  be  measured  and  used  to  predict 
the  characteristics  of  the  body  causing  the  trail,  essentially  its  shape 
and  weight. 

The  available  theoretical,  and  experimental  literature  is  reviewed  and 
the  basic  aspects  of  hypersonic  trails  are  presented.  In  the  case  of  ther¬ 
modynamic  equilibrium,  a  universal  solution  is  found  for  the  velocity  and 
enthalpy  distributions  at  a  station  behind  the  body  where  the  pressure  has 
reached  its  ambient  free -stream  value.  This  solution  is  given  in  terms  of 
the  coordinate  defining  the  shape  of  the  bow  shock  wave.  It  is  found  that 
the  nondimens ional  velbcity  and  enthalpy  profiles  depend  strongly  on  the 
drag  coefficient  alone.  These  analytical  results  are  in  good  agreement 
with  available  numerical  solutions  by  the  method  of  characteristics. 

The  thermal-conduction  part  of  the  trail  is  .Iso  studied.  An  analytic 
solution  is  found  for  the  case  of  variable  thermal  conductivity.  The  length 
of  the  trail  based  on  a  minimum  ionization  level  is  calculated  at  different 
altitudes  for  an  illustrative  re-entry.  It  is  shown  that  to  within  a  good 
approximation  this  length  is  directly  proportional  to  the  local  atmospheric 
density,  to  the  drag  coefficient,  and  to  the  object’s  cross-sectional  area, 
for  a  constant  flight  velocity  and  f6r  a  relatively  blunt  body  with  bound 
layer  effects  neglected.  The  influence  of  the  trailing  shock  on  the  cor  c- 
tion  part  of  the  trail  is  discussed.  A  preliminary  study  is  also  made  ' 
the  trail  under  chemically  frozen  conditions. 
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SYMBOLS 

Cp  ■=  drag  coefficient  (forebody  form  drag) 

Cp  =  specific  heat  at  constant  pressure 
D  =  diffusion  coefficient 
d  =  diameter 

E  a  electric  field  intensity 
f  =  velocity  ratio  defined  in  Eq.  (15) 
h  =  enthalpy 

E  =>  non  dimensional  enthalpy  defined  in  Eq.  (6) 

J,K  =  Bessel  functions 

k  =  coefficient  of  proportionality  defined  in  Eq.  (3) 
k-(.  =  thermal  conductivity 

L  =  length  of  trail 
M  =  Mach  number 

m  =  mass,  or  exponent  as  defined  in  Eq.  (3) 

N  =  concentration  of  a  species 
Pr  =  Prandtl  number 
p  =  pressure 
Re  =  Reynolds  number 
R(x)  =  function  defining  bow  shock  wave 
Ry  =  Howarth  variable  defined  in  Eq.  (l4) 
r  =  radial  distance  of  a  streamline  defined  in  Fig.  1 
S  =  entropy 
T  =  temperature 
t  =  time 


U 


velocity  in  the  axial  direction 
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X  =  distance  defined  In  Eq.  (31) 
x  =  axial  distance 

GREEK  SYMBOLS 

a(x)  =  nondlmenslonal  enthalpy  at  the  axis  defined  in  Eq.  (34) 
a*  =  numerical  coefficient  given  in  Ref.  18 
(3  =  Gaussian  depth  defined  in  Eq.  (l8) 

3'  •-=  Gaussian  depth  defined  in  Eq.  (34) 

7  =>  ratio  of  specific  heats 

71  *»  numerical  coefficient  given  in  Ref.  16 

7^  “  equivalent  ratio  of  specific  heats  defined  in  Eq.  (7) 

5  =>  effective  thickness  of  conduction  part  of  trail 

5^  =  boundary -layer  thickness 

8g  =  shock-wave  angle  (see  Fig.  l) 

6  =  characteristic  thickness  of  the  wake 

w 

5  =  electron  attachment  coefficient 

p  =  mobility  or  viscosity 

p  =  density 

cjj  =  exponent 

SUBSCRIPTS 

E  =  corresponds  to  equilibrium  flow 

F  =  corresponds  to  frozen  flow 

s  =  corresponds  to  stagnation  condition 

x  =  distance  as  given  by  Fig.  1 

o 

00  =  corresponds  to  the  free  stream 
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I,  INTRODUCTION 

Trails  of  bodies  moving  at  hypersonic  speeds  are  of  interest  to  meteor 
physicists  and  those  concerned  with  the  problems  of  describing  man-made  ob¬ 
jects  re-entering  the  earth's  dense  atmosphere. 

Several  papers  have  appeared  recently  that  attempt  to  give  an  account 
of  the  different  phenomena  involved  in  the  production  of  the  trail  and  their 
possible  use  in  predicting  the  geometry  and  flight  characteristics  of  the 
object  causing  it.  For  those  interested  in  acquiring  a  knowledge  of  the 
fundamental  behavior  of  hypersonic  trails  under  different  flow  regimes,  a 
paper  by  Feldman  offers  an  excellent  introduction.^ 

Figure  1  depicts  the  different  regions  of  the  flow  around  a  blunt  body. 
The  free-streara  region  Is  denoted  by  Aj  B  is  the  subsonic  region  behind  the 
strong  part  of  the  shock  wave;  C  corresponds  to  the  inviscid  region  between 
the  bow  shock  wave  and  the  trailing  shock.  The  boundary  layer  over  the 
body  is  denoted  by  E,  whereas  the  "near  wake"  is  shown  as  region  F.  The 
part  of  the  flow  that  is  shed  behind  the  body  by  passing  through  the  boun¬ 
dary  layer  forms  aa  inner  core,  G,  which  will  be  referred  to  here  as  the 
wake.  After  expanding  in  region  C,  the  flow  is  recomprcssed  through  the 
trailing  shock  entering  region  D  in  order  to  be  re-expanded,  thi3  time  in 

the  direction  of  flow  x.  At  a  station  x  =  x  the  pressure  will  have  reached 

o 

the  free-stream  value  and,  from  there  on,  thermal  conduction  and  other 

transport  mechanisms  will  determine  the  rate  of  the  energy  diffusion.  This 

last  region,  H,  is  normally  much  longer  than  the  distance  x  and  forms  what 

o 

we  shall  call  the  trail.  To  summarize,  viscous  and  heat-conducting  effects 
arc  prominent  in  regions  E,  F,  and  G.  Regions  B,  C,  and  D  are  considered 
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inviscid  and  adiabatic,  whereas  in  region  H  viscous  effects  are  negligible 
compared  to  conduction. 

Hie  picture  Just  given  refers  to  a  blunt  body  for  which  the  trail  is 
essentially  the  result  of  the  large  entropy  generation  created  by  the  detached 
shock  wave.  It  is  obvious  that  slender  pointed  bodies,  producing  weaker 
shocks,  will  produce  fewer  thermodynamic  observables  in  the  trail.  In  the 
blunt  geometry,  the  exact  shape  of  the  afterbody  is  not  expected  to  influence 
appreciably  the  characteristics  of  the  flow  in  the  trail.  The  validity  of 
this  statement  hinges  essentially  around  the  degree  of  irreversibility 
created  by  the  trailing  shock;  inasmuch  as  sharply  closing  streamlines  are 
not  expected  to  form  behind  the  body,  the  discontinuity  due  to  the  second 
shock  may  be  assumed  reversible  on  a  first  approximation.  This  point  will 
be  elaborated  later. 

As  a  first  rough  approximation,  the  trail  will  depend  on  the  following 
parameters:  altitude,  velocity,  and  a  characteristic  length.  The  following 
a  priori  qualitative  observations  regarding  the  nature  of  the  trail  are 
based  on  general  aerodynamic  and  thermodynamic  arguments.  At  low  altitudes, 
where  the  density  is  fairly  high,  we  expect  the  flow  to  be  in  thermodynamic 
equilibrium  everywhere.  We  also  expect  that  the  boundary -layer  thickness 
around  the  body  will  be  small  and  that,  as  a  result,  a  small  percentage 
of  the  energy  associoted  -with  the  flow  will  go  through  the  boundary  layer. 

In  consequence,  the  viscous  wake  G  is  expected  to  be  thin  with  a  low 
energy  content  compared  to  the  total  energy.  In  such  a  case,  we  may  neglect 
the  existence  of  regions  E,  F,  and  G  (see  Fig.  l). 

Consider  now  the  possibility  of  turbulence.  It  seems  that  low 
critical  Reynolds  numbers  are  needed  in  order  for  the  viscous  wake  G  to 
become  unstable.  The  influence  of  turbulence  on  the  observables  in  the 
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trail  will  be  more  important  at  those  altitudes  where  more  energy  will  be 
associated  with  a  thicker  boundary  layer.  On  the  other  hand,  the  effect  of 
region  G,  even  though  turbulent,  will  be  minimized  if  its  rate  of  growth 
in  the  radial  direction  is  small.  We  conclude  that  at  low  altitudes,  with 
a  thin  boundary  layer,  only  the  shape  of  the  induced  detached  shock  will 
determine  the  nature  of  the  trail,  and  that  knowledge  of  this  shape  should 
be  adequate  for  the  description  of  the  trail. 

At  high  altitudes,  where  the  density  is  low,  relaxation  effects  will 
be  important.  In  the  subsonic  region  B,  where  the  densities  are  high  due 
to  the  normal  character  of  the  shock,  the  flow  will  be  dissociated,  but  in 
thermodynamic  equilibrium.  As  the  flow  expands  to  lower  densitites  the 
process  of  recombination  of  the  dissociated  species  lags  behind,  since  it 
is  governed  by  a  three-body  collision;  nonequilibrium  effects  then  become 
important.  In  this  case,  the  degree  of  departure  from  equilibrium  will  be 
governed  net  only  by  the  local  state  of  the  fluid,  but  also  by  a  characteristic 
length.  For  the  same  velocity,  smaller  lengths  will  be  traveled  at  shorter 
times  and,  hence,  a  smaller  number  of  collisions  will  be  associated  with 
the  motion  of  the  particles,  resulting  in  larger  departures  from  equilibrium. 
The  phenomenon  becomes  more  complicated  when  the  rate  of  recombination  of 
some  species  is  different  from  that  of  others  (for  instance,  atom-atom 
recombination  versus  ion-electron). 

In  the  final  analysis,  refinement  of  the  calculations  will  depend  on 
the  nature  of  the  observable.  For  instance,  if  observation  of  a  trail  in 
thermodynamic  equilibrium  is  planned  with  a  radar  team,  the  electron- 
concentration  profiles  are  of  importance;  some  error  in  the  prediction  of 
the  thermodynamic  state  (essentially  temperature)  can  be  tolerated. 
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On  the  other  hand,  if  the  radiant  intensity  of  the  trail  is  of  interest, 
then  because  of  the  stronger  dependence  of  radiation  on  temperature,  a 
small  error  in  temperature  may  lead  to  a  manifold  error. In  this 
connection  it  should  be  emphasized  that  the  definition  of  the  trail  length 
depends  on  the  nature  of  the  observable  and  on  the  observing  instrumen¬ 
tation.  For  a  chosen  observable,  the  trail  length  may  be  defined  as  the 
distance  behind  the  body  at  which  a  given  instrument  is  unable  to  dis¬ 
criminate  or  record  from  ambient  conditions. 

Chemical  relaxation  effects  become  important  above  an  altitude  of 
about  100,000  ft  for  objects  of  reasonable  size  re-entering  at  satellite 
or  ICBM  velocity.  At  an  altitude  roughly  above  300,000  ft,  most  of  the 
diatomic  particles  are  absent,  and  me an- free -path  effects  become  important. 
Coupling  effects  in  the  ionosphere  with  the  earth's  magnetic  field  are 
possible,  but  this  region  has  yet  to  be  studied  in  a  systematic  way.  It 
is  obvious,  nevertheless,  that  the  observables  are  no  longer  associated 
with  the  existence  of  an  entropy  trail  in  the  sense  described  in  this 
section. 

A  study  of  the  hypersonic  trail  aims  at  the  development  of  thermo¬ 
dynamic  profiles  of  those  parameters  which,  with  our  present  state  of 
knowledge,  can  be  measured.  These  measurements  are  used  in  an  attempt  to 
predict  the  characteristics  of  the  body  causing  the  trail,  essentially 
its  shape  and  weight. 
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II.  DISCUSSION  OF  AVAILABLE  RESULTS 

THE  TRAIL  IN  THERMQDY NAMI C  EQUILIBRIUM 

For  the  reasons  enumerated  in  the  previous  section,  low  altitudes  will 
be  associated  with  thermodynamic  equilibrium;  furthermore,  the  role  of 
viscosity,  confined  to  regions  E,  F,  and  G  (see  Fig.  l)  will  be  relatively 
small. 

Assuming  that  the  amount  of  fluid  entering  the  region  G  is  equal  to 
the  fluid  mass  going  through  the  boundary  layer  over  the  body,  we  have  the 
following  order -of -magnitude  equality: 

Here  is  the  characteristic  thickness  of  the  woke  G  and  5^  the  boundary- 
layer  thickness  over  the  body.  We  conclude  that 

6  Jr 

v  o  v  b'  c 

This  means  that  as  long  as  6^  remains  small,  5  will  also  be  small.  The 

°  b  ’  w 

same  argument  was  offered  in  Ref.  1. 

The  calculation  of  the  trail  proceeds  in  the  following  way.  For  given 

geometry  and  flight  conditions,  we  calculate  the  shape  of  the  shock  wave 

and  the  inviscid  flow  behind  it  by  making  use  of  the  method  of  characteristics. 

We  will  start  with  a  blunt  nose  and  consider  different  afterbodies  and 

(l  2) 

trailing  conditions.  Feldman  '  has  investigated  the  case  of  a  two- 
dimensional  blunt  nose  followed  by  a  flare  of  an  angle  between  0  and  30  deg 
and  by  an  expansion  in  the  near-wake  region  F  between  0  and  20  deg.  These 
calculations  were  carried  out  for  both  a  short  and  a  long  body.  Similar 
calculations  have  been  carried  out  for  a  hemisphere  followed  by  a  frustum, 
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or  a  true  streamline  and  then  a  cylinder.*  Frustum  or  free-streamline 
angles  up  to  30  deg  were  considered.  All  of  the  above  results  indicate 
that  at  a  distance  of  about  10  radii  behind  the  body,  the  thermodynamic 
state  is  very  nearly  the  same  as  the  one  corresponding  to  an  isentropic 
flow  for  expansion  angles  up  to  20  deg.  A  difference  in  temperature  of 
about  10  per  cent  was  found  for  a  frustum  of  30  deg,  but  it  should  be 
recognized  that  this  is  an  extreme  case  not  to  be  anticipated  in  reality. 
An  expansion  angle  of  20  deg  is  supported  by  available  experimental 
evidence.**  The  effect  of  the  exact  geometry  of  the  afterbody  on  the 
equilibrium  electron  concentration  and  radiation  emitted  in  the  optical 
(  ^  2,000  A°  to  10,000  A°)  and  infrared  (  =  5-3  microns)  region  of  the 
spectrum  has  been  evaluated  in  Ref.  1.  The  final  conclusion  is 

—  [the]  magnitude  of  thermal  radiation  and 
electron  density  that  are  arrived  at  on  the  basis  of 
calculations  of  a  simple  blunt  body  geometry  like  a 
hemisphere— cylinder,  that  does  not  include  the  geo¬ 
metrical  details  of  the  afterbody,  may  be  off  by  a 
factor  of  not  more  than  3- 

Since  the  reflecting  qualities  of  a  radar  beam  over  an  ionized  trail 
depend  on  its  plasma  frequency,  which  is  directly  proportional  to  the 
square  root  of  the  electron  concentration,  the  error  in  the  plasma 


frequency  is  diminished.  In  this  memorandum,  the  problem  of  interaction 

between  a  radar  beam  generating  a  field  E  varying  as  exp(ix  x  +  ico^t)  and 

,  7 

a  medium  containing  II  electrons/crcr'  will  not  be  covered.  Hie  following 


Personal  communication  with  L.  Hromas,  Space  Technology  Laboratories, 

Inc. 

v  This  expansion  angle  is  eventually  determined  by  consideration  of 
viscous  effects.  No  theoretical  work  seems  to  be  available  at  present  in 
this  region. 


RM-2682 -1 
9 


For  rOp  >  the  signal  cn^  is  reflected,  whereas  for  cu^  <  av^  the  radar 
wave  is  damped  out.  The  electron  densities  occurring  in  the  ionosphere 
are  of  the  order  of  10^  electrons/cmA  It  follows  that  no  trail  with  a 
lower  electron  concentration  can  be  detected  by  a  radar  beam. 

Having  established  the  reversible  character  of  the  geometry 
behind  the  blunt  nose,  the  inviscid  calculation  is  performed  over  a 
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f  i  o ) 

hemisphere-cylinder  configuration.  Feldman v  ,l~  gives  all  the  thermo¬ 
dynamic  profiles  for  a  hemisphere-cylinder  geometry  extending  to  the 
station  x  =  xq  of  Fig.  1  where  the  pressure  everywhere  has  reached 
approximately  its  free-stream  value.  For  an  altitude  of  60,000  ft  and 
velocity  of  17,500  ft/ sec  the  pressure  reaches  its  free-stream  value  at 
a  distance  of  about  5 6  radii. For  an  altitude  of  100,000  ft  and 
velocities  of  15,000,  20,000,  and  25,000  ft/sec,  the  distances  are  about 
40,  80,  and  100  radii,  respectively.*  In  all  cases,  this  is  a  small 
fraction  of  the  total  length  of  the  trail  over  which  temperature  and 
electron-concentration  gradients  will  persist. 

As  was  stated  in  the  previous  section,  from  the  station  x  =  xq  onwards, 
the  trail  will  cool  off  after  tennination  of  the  adiabatic  expansion  by 
the  mechanism  of  heat  conduction.  Viscous  dissipation  offers  a  very  small 
contribution,  since  the  velocity  deficiency  at  the  station  x  =  xq  at  the 
center  line  is  of  the  order  of  20  per  cent  of  the  free-stream  value. 
Solutions  for  the  conduction-controlled  part  of  the  trail  are  available 
in  Refs.  1-6. 

In  Ref.  1,  the  initial  enthalpy  and  velocity  profiles  found  by  the 

method  of  characteristics  and  corresponding  to  the  station  x  =  x  are 

2,  ° 

approximated  by  a  Gaussian  distribution  of  the  form  e  where  R^  is 

related  to  the  physical  radius  r  through  the  Hovarth  transformation 

4  =  2Jf~  rdr 

o 


*  Personal  communication  with  S.  Feldman,  Electro-Optical  Systems,  Inc. 
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For  this  purpose,  since  all  the  characteristic  profiles  are  calculated  at 

the  surface  of  the  cylinder,  the  value  at  r  =  rQ  was  transferred  to  r  =  0 

and  the  interval  between  the  center  line  and  r  =  r  was  faired  with  a 

o 

smooth  curve.  For  the  evolution  of  the  trail  during  its  thermal  mixing, 
an  integral  method  was  applied,  satisfying  the  exact  boundary  conditions 
and  differential  equation  at  r-»0,  and  also  the  condition  that  the  total 
energy  and  momentum  in  the  radial  direction  extended  to  infinity  remain 
constant  at  all  stations  x.  From  the  final  solution,  it  was  found  that 
the  error  made  by  neglecting  the  transfer  of  momentum  in  detail  was  very 
small.  After  the  thermodynamic  state  inside  the  trail  was  thus  established, 
calculations  were  made  of  electron-concentration  profiles.  The  decay  of 
infrared  and  optical  radiation  was  also  calculated  along  the  axial  stream¬ 
line.  The  results  are  given  for  altitudes  of  100,000  and  250,000  ft  with 
velocities  ranging  between  15,000  and  35>000  ft/sec. 

In  Ref.  3  a  solution  of  the  conduction-controlled  trail  is  offered, 
the  equations  of  conservation  being  solved  by  the  method  of  finite 
differences  with  the  help  of  a  digital  computer.  A  comparison  with 
Feldman's  integral  method^^  showed  that  the  trail  lengths  based  on 
electron-concentration  profiles  were  about  half  the  length  of  those 
obtained  by  the  finite-difference  method.  This  discrepancy  has  been 
attributed  in  Ref.  3  partly  to  the  fact  that  the  Gaussian  initial  profiles 
chosen  by  Feldman  did  not  fit  away  from  the  axis  when  compared  with  those 
he  obtained  by  the  method  of  characteristics  (see  Fig.  2).  Actually, 
comparison  of  the  initial  profile  chosen  by  the  Goulards  in  the  physical 
plane^"^  with  the  ones  obtained  by  Feldman^  shows  that  the  discrepancy 
is  due  to  a  presumably  erroneous  interpolation  in  Ref.  3  between  the 


Comparison  of  initial  enthalpy  profiles  before  thermal  mixing 
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velocities  of  20,000  and  ?5,000  ft/ sec.  Inspection  of  Fig.  2  shows  that 
Feldman's  predictions  of  the  trail  length  should  be  closer  to  the  accurate 
ones.  Furthermore,  the  findings  reported  in  Section  III  will  theoretically 
justify  Feldman's  choice  of  the  Gaussian  widths  for  his  solution. 

An  equilibrium  solution  for  a  blunt  body  for  an  entropy-induced 
shock,  wave  is  also  presented  in  Ref.  4.  The  initial  velocity  profile 
corresponding  to  the  section  x  =  xq  is  calculated  by  the  streamtube 
method,  starting  with  a  known  shock-wave  shape  and  using  real-gas  tables 
for  the  isentropic  expansion  behind  the  body.  The  velocity  profile  found 
is  highly  irregular,  and  the  authors  approximate  it  by  using  the  method  of 
least  squares  with  a  third-degree  polynomial,  exhibiting  the  same  behavior 
as  a  Gaussian  distribution.  The  deviations  from  this  approximation  are 
large  but  no  explanation  of  the  reason  is  offered  in  the  paper.  A 
momentum-integral  method  is  then  used  for  the  study  of  the  trail.  The 
enthalpy  is  calculated  from  the  Crocco  integral,  the  Prandtl  and  Lewis 
numbers  both  being  assumed  equal  to  unity.  The  electron  concentration  is 
given  along  the  center  line  for  an  altitude  of  200,000  ft  and  a  velocity 
of  23,000  ft/sec. 

RELAXATION  EFFECTS 

At  present  there  is  no  definitive  quantitative  analysis  available  in 
the  literature  on  nonequilibrium  effects,  apart  from  some  preliminary 
calculations . 

In  Ref.  7  a  study  is  made  of  the  importance  of  relaxation  effects 
along  the  stagnation  streamline  over  a  blunt  configuration.  Flight 
velocities  of  15,00u  and  25,000  ft/sec  with  altitudes  between  154,000  and 
246,000  ft  have  been  considered.  The  major  conclusion  is  that  for  altitudes 
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higher  than  150,000  ft  the  flow  in  the  inviscid  supersonic  expansion  region 
of  a  blunt  body  will  be  essentially  chemically  frozen.  It  should  be  kept  in 

/  a\ 

mind  that  the  transition  to  frozen  conditions  progresses  very  rapidly. 

A  thorough  calculation  of  the  whole  inviscid  field  should,  of  course, 
include  a  study  of  all  the  streamlines  crossing  the  bow  shock  wave. 

In  Ref.  3,  it  was  found  that  transition  to  frozen  flow  develops  at  about 
45  deg  aft  of  the  stagnation  point  for  a  velocity  of  23,000  ft/sec,  an  alti¬ 
tude  of  100,000  ft,  and  a  characteristic  length  of  1  ft.  It  was  also  calcu¬ 
lated  that  for  the  same  flight  velocity  and  same  body  diameter  the  transition 
point  moves  closer  to  the  stagnation  region  at  higher  altitudes.  For  a 
velocity  of  23,000  ft/sec  and  a  radius  of  1  ft,  transition  is  almost  at  the 
stagnation  point  for  altitudes  above  150,000  ft.  This  means  that  above 
these  altitudes  the  calculation  could  be  simplified  by  assuming  frozen 
conditions  almost  immediately  behind  the  shock  wave.  It  should  be  added 
that  this  procedure  will  be  more  correct  for  the  streamlines  closer  to  the 
stagnation  region  where  most  of  the  irreversibility  (and  interest)  lies. 

As  a  first  approximation  for  the  calculation  of  the  temperature  profile 
at  the  region  x  =  xq,  on  the  assumption  of  complete  frozen  conditions,  one 
could  use  the  equivalent  effective-specific-heat  ratio  for  frozen  air  as 
given  in  Ref.  9  or  10.  In  general,  this  effective  value  approaches  the 
ratio  of  5/3  of  a  monatomic  gas.  The  temperature  at  the  axial  streamline 

at  the  station  x  =  x  would  be  lower  if  transition  to  frozen  conditions 

o 

occurred  earlier.  This  is  self-evident  since  the  isentropic  expansion  at 

a  higher  value  of  7„  takes  place  over  a  greater  length  of  time.  On  the 
r 

other  hand,  we  expect  to  find,  for  the  same  body  diameter  and  altitude, 
smaller  dependence  of  the  temperature  at  the  station  x  =  xq  on  the  flight 
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velocity,  as  transition  moves  to  the  stagnation  point.  This  may  be  under¬ 
stood  as  follows:  The  value  of  7  is  higher  for  higher  compressibility 

r 

factors  which,  in  turn,  are  higher  for  higher  flight  velocities  or  flight 
Mach  numbers  Since 


T  ~ 
x 

o 


T 

s 


7 y  ^-/7 y 


the  increases  in  M  will  be  moderated  by  decreases  in  the  ratio  2/V^. 

CO  r 

These  qualitative  remarks  are  substantiated  by  the  calculations  in  Ref.  1, 
undertaken  for  an  altitude  of  200,000  ft  and  flight  velocities  between 
10,000  and  30,000  ft/sec  (see  Fig.  7  of  Ref.  1). 

The  variation  of  the  electron  concentration  behind  the  body  will  next 
be  considered.  The  rate  of  recombination  of  the  ions  with  electrons  is  of 
primary  importance  here.  It  is  evident  that  sxower  ion-electron  recombina¬ 
tion  (high  altitudes)  will  result  in  an  electron  concentration  behind  the 
body  that  i6  higher  than  the  concentrations  obtained  for  thermodynamic 
equilibrium.  Nevertheless,  at  intermediate  altitudes,  with  slow  atom-atom 
recombination  and  fast  ionic  recombination,  it  is  conceivable  that  lower- 
than-equilibrium  electron  concentrations  may  result,  since  under  these 
conditions  the  electron  density  will  be  in  equilibrium  with  the  local 
frozen  temperature  (which  is  lower  than  the  equilibrium  temperature). 

This  was  pointed  out  first  by  S.  C.  Lin.'1'' 

Frozen-flow  calculations  of  electron  concentrations  along  the  axis  of 
the  trail  are  elso  reported  in  Ref.  U.  Again,  as  in  the  equilibrium 
calculations,  all  the  profiles  at  x  =  xq  computed  by  assuming  frozen 
conditions  immediately  behind  the  shock  wave  were  highly  irregular.  For 
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fitting,  third-degree  polynomials  were  used.  The  study  consists 
essentially  in  a  conservation  and  redistribution  of  the  electron  con¬ 
centration  in  the  trail.  The  exact  thermochemistry  of  the  phenomenon  is 
not  taken  into  account. 

A  summary  of  the  work  of  Bortner  is  given  in  Ref.  5*  Dissociation 
and  ionization  reactions  have  been  taken  into  account.  It  is  reported 
that  an  order-of-magnitude  change  in  the  oxygen-dissociation-rate  constant 
can  cause  increases  by  factors  of  over  100  in  the  electron  density.  A 
10-deg  spherically  blunted  cone  at  an  altitude  of  160,000  ft  is  considered. 

The  study  is  made  along  the  stagnation  streamline.  It  is  found  that  the 
electron  concentration  six  characteristic  lengths  behind  the  body  is  one 
order  of  magnitude  higher  if  the  flow  is  frozen.  On  the  other  hand,  lower- 
than-equilibrium  electron  concentrations  are  predicted  for  mixed  nonequilibrium 
flow.  (This  is  the  trend  anticipated  in  Ref.  1  and  discussed  above.)  A 
study  is  also  made  of  contamination  effects  in  the  electron  concentration. 
Equilibrium  is  assumed  and  contamination  of  one  part  sodium  atoms  per 
million  and  one  part  sodium  atoms  per  ten  thousand  parts  of  equilibrium 
air  is  considered.  The  electron  concentration  is  increased  by  at  least 
three  orders  of  magnitude  for  the  altitude  of  250,000  ft  and  Mach  number  2k. 

All  of  the  above  calculations  extend  to  a  few  characteristic  lengths  behind 
the  body  in  the  inviscid  region  and  do  not  cover  the  trail  proper  (region  H 
of  Fig.  1). 

At  rather  high  altitudes,  where  low  densities  prevail,  the  mechanism 
of  ambipolar  diffusion  is  of  importance.  Tnis  mechanism  occurs  because 
the  rate  of  diffusion  of  the  ions  and  electrons  in  the  gas  is  different, 
thus  causing  a  space-charge  field.  Consideration  of  the  particular  law  of 
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motion  for  every  species  is  difficult;  nevertheless,  a  "gross"  approach  is 
traditionally  used.  We  assume  that  the  ion  and  electron  concentrations 
remain  the  same  and  that  their  average  velocities  are  equal.  With  the 
subscripts  i  and  e  corresponding  to  ions  and  electrons,  and  under  the 
assumption  that  the  gross  space  charge  is  represented  by  an  electric  field 
of  magnitude  E,  we  write 


N  v 
i  i 


+  n1  E 


N  v 
e  e 


=  -D 

e 


dU 

e 

dx 


E 


This  is  a  charge-conservation  equation  for  the  motion  in  direction  x.  The 
first  term  on  the  right  represents  the  diffusion  contribution,  whereas 
the  second  one  gives  the  contribution  from  a  velocity  imparted  to  the 
particle  by  the  presence  of  the  space  charge  represented  by  E.  This 
velocity  is  proportional  to  the  mobility.  Eliminating  E  from  the  above 
two  equations  by  making  Ng  =  =  N,  v^  =  vg  =  v,  we  have 


N  v  = 


Vi  +  Ve\ 

^e  +  ^e  / 


dN 

dx 


which  defines  the  "coefficient  of  ambipolar  diffusion"  D&  as 


D 

a 


Ve  *  Ve 

“i  *  “e 


Simple  kinetic  theory  shows  that,  irrespective  of  the  difference  of  mass, 

D  /p  kT/e,  where  k  is  Boltzmann's  constant;  taking  into  account 

the  fact  that  p  »  p .  ,  we  find 
e  l 
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Another  mechanism,  important  under  a  wide  range  of  conditions,  is  the 
attachment  of  electrons  to  neutral  species,  thus  creating  negative  ions; 
in  the  case  of  air,  this  would  involve  the  attachment  of  electrons  to 
diatomic  oxygen.  The  rate  of  attachment  is  proportional  to  the  density  of 
the  electron  and  to  the  square  of  air-molecule  density. 

The  electric-charge-continuity  equation  takes  the  form 

-  SjJ  «  CNN2  +  a  N2  -  D 
ot  a  r  a  v 

The  first  term  on  the  right  is  the  contribution  through  attachment,  the 
second  through  recombination,  and  the  third  through  ambipolar  diffusion; 

N&  is  the  air-molecule  density. 

The  mechanisms  of  ambipolar  diffusion,  recombination,  and  attachment 
are  considered  in  some  detail  in  Ref.  12.  The  computations  are  performed 
for  a  5-da.  re-entry  body.  Temperature  variations  in  the  reaction-rate 
coefficients  have  been  neglected,  since  the  variations  in  density  with 
altitude  are  much  more  important.  It  is  found  that  recombination  prevails 
over  ambipolar  diffusion  for  all  altitudes  under  260,000  ft,  above  which 
diffusion  is  the  predominant  mechanism.  On  the  other  hand,  attachment  is 
found  to  be  considerable  in  all  cases  if  the  high  values fbr  £  measurement 
in  drift  tubes  are  adopted.  It  is  pointed  out  that  during  an  actual 
re-entry  much  lower  attachment  rates  will  be  observed;  this  is  because 
drift  measurements  disturb  the  gas  very  little,  whereas  during  actual 
re-entry  the  air  is  disturbed  considerably.  Electron  profiles  for 


characteristic  lengths  of  1  cm,  5  in.,  and  5  ft  at  a  velocity  of  25,500 


ft/sec  have  been  calculated  for  a  re-entry  example. 


Comparisons  are 


made  with  the  attachment  mechanism  suppressed.  The  absence  of  rhis 


mechanism  alters  significantly  the  electron  concentrations,  especially 
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at  the  lower  altitudes  where  the  dependence  of  is  felt  more  strongly 
in  the  diffusion  equation. 


THE  PRESENCE  OF  TURBULENCE 

In  Ref.  14  the  problem  of  the  hydrodynamic  stability  of  a  Gaussian 
velocity  profile  in  an  incompressible  viscous  wake  is  analyzed  by  using 
the  method  of  Galerkin  for  the  solution  of  the  Sommerfeld— Orr  equation. 

It  is  found  that  the  critical  Reynolds  number  for  antisymmetric  perturbations 
begins  to  develop  at  the  value  19  with  the  Reynolds  number  defined  from 
the  relation 


iReT 


U  dC_ 
oo  D 


crlt 


v  2  n 

where  d  is  defined  as  the  effective  cross-sectional  dimension  of  the  body 
causing  the  viscous  wake.*  The  experiments  of  Ref.  15,  undertaken  for 

( i>5  \ 

wakes  behind  cylinders,  confirm  the  findings  of  the  theory.  Kovasznayx  •" 
made  his  measurements  in  the  low-Reynolds— number  range  and  observed  that 
the  familiar  KArraAn  vortex  streets  appear  only  above  Re  =  40  and  that  they 
remain  stable  and  regular  for  Reynolds  numbers  below  160.  In  Roshko's 
paper, ^  ^  wakes  were  observed  for  Reynolds  numbers  between  40  and  10,000. 
His  experimental  results  coincide  with  Kovasznay's  for  Reynolds  numbers 
between  40  and  150  where  he  observed  that  the  free  vortices  move  downstream, 
decaying  by  "the  laminar  viscous  mechanism.  For  Reynolds  numbers  greater 
than  300  the  vortices  diffused  downstream  turbulently  and  the  wake  became 
fully  turbulent  in  about  40  to  50  diameters  downstream.  The  region  of 


^Symmetric  perturbations  exhibit  instability  at  (Re)„r^.j.  55;  the 

regular  Reynolds  numbers  are  about  twice  as  high. 
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Reynolds  numbers  between  150  and  300  was  a  transition  range  in  which 
irregular  velocity  fluctuations  occurred. 

No  experiments  are  available  in  the  case  of  three-dimensional  bluff 
bodies  at  hypersonic  speeds,  hence  the  extent  of  the  influence  of  three- 
dimensional  and  compressibility  effects  in  the  development  of  turbulence 
inside  the  viscous  wake  remains  unknown. 

On  the  other  hand,  from  elementary  wake  theory  it  is  known,  that  the 
effective  width  of  the  wake  increases  in  the  direction  of  flow  x  with  a 
power  between  l/2  and  l/ 3  for  two-dimensional  and  circular  wakes,  respec¬ 
tively.  This  means  that  if  the  wake  contains  a  small  amount  of  the  total 
energy,  as  is  the  case  for  thin  boundary  layers  (low  altitudes),  the 
spreading  of  the  turbulent  regime  towards  the  outside  in  the  laminar 
region  E  of  Fig.  1  will  not  cause  significant  change  in  the  trail  lengths 
computed  by  considering  laminar  flow.  Conversely,  at  higher  altitudes 
where  a  considerable  portion  of  the  energy  associated  with  the  moving  body 
will  go  through  a  thick  boundary  layer,  turbulence  will  spread  rather 
rapidly,  and  the  results  of  any  laminar  calculation  will  give  a  trail 
length  much  longer  than  the  actual  one.  It  should  be  emphasized  that  there 
is  no  unique  definition  of  the  trail  length.  Depending  on  the  chosen 
observable,  it  may  be  defined  as  the  distance  behind  the  body  at  which  a 
given  instrument  is  unable  to  discriminate  or  record  from  ambient  conditions. 

In  Ref.  6  a  calculation  is  performed  for  the  turbulent  trail.  The 
velocity  and  enthalpy  profiles  corresponding  to  station  x  =  xq  were 
estimated  by  assuming  an  elliptical  shock  wave  and  an  isentropic  expansion 
to  atmospheric  pressure  along  each  streamline.  The  semiempirical  eddy- 
diffusivity  theory  for  the  calculation  of  the  turbulent  shear  was  used, 
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(17 ) 

based  on  the  work  of  Townsend.  It  was  assumed,  furthermore,  that  the 

entire  trail  becomes  turbulent  immediately  after  the  station  x  =  the 

trail  is  therefore  found  to  be  much  shorter  than  the  laminar  one  as  computed 
by  Feldman^ ^  for  the  same  flight  conditions  (60,000-ft  altitude  and 
17,500— ft/sec  velocity).  For  the  reasons  explained  previously,  the  assumption 
of  turbulence  in  every  streamline  does  not  seem  to  be  realistic,  especially 
for  the  low  altitude  chosen  in  the  numerical  example.* 

Equilibrium  turbulent  calculations  have  also  been  performed  by  Long 
and  reported  in  Ref.  5*  A  wide  range  of  conditions  is  covered  with 
altitudes  between  150,000  and  250,000  ft  and  velocities  ranging  between 
15,000  and  25,000  ft/sec.  The  geometry  chosen  was  a  spherically  blunted 
10-deg  cone.  A  turbulent  Prandtl  number  of  unity  was  assumed.  It  appears 
that,  as  in  Ref.  6,  turbulence  has  been  assumed  to  prevail  both  in  the 
viscous  wake  and  the  trail.  Boundary -layer  effects,  however,  have  been 
neglected. 

*In  a  later  publication  W.  W.  Short  and  A.  Hoehstim  have  retracted  this 
assumption.  See  The  Effect  of  Additives  in  the  Turbulent  Wake  of  a  Re-Entry 
Vehicle ,  Convair,  Physics  Section  Ph-0c’7 -M,  August  1000. 
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III.  APPROXIMATE  UNIVERSAL  SOLUTIONS 


The  solutions  presented  in  the  literature  and  discussed  in  the  previous 
section  have  been  obtained  partially  through  numerical  calculations  performed 
by  digital  computers.  Given  a  certain  geometry  and  flight  conditions,  a 
numerical  answer  is  provided  on  the  nature  of  the  trail  without  an  intimate 
knowledge  of  the  basic  coupling  mechanism  of  the  different  parameters 
involved.  These  parameters  remain  buried  inside  the  numerical  computations, 
and  it  is  impossible  to  retrieve  them  in  a  clear  form  from  the  final  result. 

In  this  section  an  attempt  is  made  to  present  the  problem  of  the 
hypersonic  trail  from  a  fundamental  point  of  view  by  using  expressions  in 
closed  form,  where  their  physical  significance  is  more  apparent.  It  will 
be  shown  that  the  closed-form  results  are  very  close  to  the  ones  predicted 
by  machines. 

THE  EXPANSION  PART  OF  THE  TRAIL 

We  first  consider  the  case  of  hypersonic  trails  in  thermodynamic 
equilibrium.  We  shall  assume,  in  the  light  of  the  discussion  in  the 
previous  sections,  that  no  other  irreversibility  is  present  in  the  flow, 
except  for  the  detached  shock  wave  ahead  of  the  body.  Let  the  shape  of 
the  shock  wave  be  given  by  an  expression  of  th?  form  R  =  R(x),  the  function 
R  being  continuous  and  admitting  a  first-order  derivative.  We  are  interested 
in  calculating  the  ratio  h/h  the  station  x  =  xq  where  the  pressure 

at  al  1  values  of  r  (see  Fig.  l)  and  all  x  >  x  is  equal  to  the  ambient 

o 

pressure  p  .  E’rom  Ref.  18  we  can  use  the  isentropic  exponents  related  to 
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real-gas  effects  in  the  temperature  range  of  3000  to  8000°K.*  The  follow¬ 
ing  sequence  of  equalities  is  self-explanatory. 


sin26) 


(y*  -l)/a*7 ' 


1  +  ^  M2  sin2  6 

2  oo 


(1) 


The  coefficients  a*  and  y '  are  given  in  Ref.  l8**  once  the  value  of  the 
entropy  immediately  behind  the  shock  wave  has  been  established.  The  value 
of  an  effective  7  for  the  Jump  across  the  shock  wave  can  be  estimated  by 
using  the  tables  of  Ref.  11  or  Ref.  19*  The  parameter  5  is  the  shock  angle 
as  indicated  in  Fig.  1,  and  can  easily  be  calculated  by  the  given  equation 
for  the  shock  wave 

sin2  6  =  - - - rr-  =  - - - *  (2) 

1  +  1/tan  6  1  +  i/(dR/dx) 

Up  to  this  point,  we  have  not  used  the  hypersonic  approximations  valid 


*From  a  Mollier  chart  we  can  easily  see  that  this  is  the  range  of 
temperatures  behind  the  shock  wave  for  all  flight  conditions  of  interest  in 
re-entry. 

**In  Ref.  18  the  symbol  a  instead  of  O’1  is  used.  In  this  memorandum 
we  reserve  a  for  later  use. 
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for  strong  shock  waves.  Equation  (l)  gives,  in  closed  form,  the  initial 
enthalpy  profile  at  the  station  x  =  xq  once  the  shape  of  the  shock  wave, 
the  altitude,  and  flight  velocity  are  given.  Let  us  consider  the  class  of 
shock  waves  given  by  power  laws  of  the  form* 

R  =  kxm  (5) 

It  is  well  known  that  such  dependence  suggests  itself  from  a  generalization 

(20  21 ) 

of  blast-wave  theories.  ’  '  (Experimental  data  also  yield  relations  of 

this  kind. ^ )  Equation  (2)  becomes 


In  the  hypersonic  approximation  and  for  all  values  of  R  smaller  than  the 
value  for  which  the  normal  Mach  number  sinB  is  strong  enough  to  Justify 
the  omission  of  unity  from  the  numerator  and  the  quantity  (7— l)/(7+l)  from 
the  denominator  of  Eq.  (l),  we  have 


*A  body  radius  of  one  unit  is  assumed  in  these  calculations  throughout. 
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In  general 


where 


(6) 


(T) 


Between  the  sane  range  of  entropies,  corresponding  to  the  range  between 
3000  and  8000°K,  the  coefficient  a'  varies  by  about  20  per  cent,  and 
the  coefficient  7'  by  about  3  per  cent,  whereas  7^  varies  by  only  3.3 
per  cent.  In  other  words,  7  is  independent  of  flight  conditions. 
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Let  us  now  examine  in  more  detail  the  case  of  a  blunt  nose  followed  by 
a  cylinder,  for  which  it  is  well  known  that  the  shape  of  the  shock  is  given 
by  a  parabola 


R  = 


(B) 


Substitution  of  Eq.  (8)  into  Eq.  (6)  gives  (for  m  =  l/2) 


h  = 


(9) 


* 


It  is  well  known  from  blast-wave  theory  that  the  coefficient  k  is 
directly  proportional  to  the  fourth  root  of  the  drag  coefficient.^^'  From 
experimental  and  theoretical  calculations  over  a  wide  range  of  configurations, 
with  drag  coefficients  varying  between  0.037  and  1.37*  Ref.  23  suggests  the 
following  correlation  for  hypersonic  flow  in  air 


R 

d 


(10) 


where,  it  will  be  remembered,  represents  the  nose  form  drag  (i.e.,  not 
base  drag  and  skin  friction). 

The  experimental  points  for  air  from  the  above  reference  are  shown  in 
Fig.  3,  along  with  a  sketch  of  the  different  configurations  tested.  In 
the  present  analysis  we  prefer  to  use  the  correlation 


R 

d 


0.90 


(11) 


Besides  the  fact  that  the  square  root  of  x  is  suggested  by  the  theory, 
the  above  relation  does  not  lack  accuracy  when  compared  with  Eq.  (10).  When 
R  and  x  are  referred  to  the  characteristic  radius  r  ,  the  coefficient  k 


*Note  the  strong  dependence  of  h  on  the  numerical  value  of  k. 


-Experimental  data  of  hypersonic  shock  shapes  correlated 
in  accordance  with  the  blast  wave  therory 
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takes  the  value 


k  =  1.27  (12) 

It  Is  worth  mentioning  that  for  an  altitude  of  100,000  ft  Feldman 
found,  through  the  method  of  characteristics,  a  shock  shape  corresponding 
to  a  value  k  =  1.25  for  all  flight  velocities  considered.  On  the  other 
hand,  for  the  altitude  of  60,000  ft  and  a  velocity  of  17,500  ft/sec,  he 
found  k  =  1.4. 

Equation  (9)  is  truly  an  expression  of  a  universal  enthalpy  profile, 
depending  solely  on  the  drag  coefficient  and  the  ratio  y  ,  both  being 

Li 

independent  of  altitude  and  flight  velocity  within  the  frame  of  hypersonic 
approximations.  It  might  be  of  interest  to  note  that  the  profile  of  Eq. 
(9)  is  analytically  the  same  as  the  one  emerging  from  a  circular  laminar 
Jet. 

The  mass-continuity  equation  at  the  two  stations  x  =  0  and  x  =»  xq 
(see  Fig.  1)  may  be  written  as 


p  U  RdR 
oo  oo 


p  U  rdr 

X  X 


o  o 


(13) 


It  is  easy  to  establish  that  the  Howarth  transformation  variable  Ra  used  in 

n 

Ref.  1  for  the  correlation  of  the  enthalpy  profiles  obtained  through  the 
method  of  characteristics  is  connected  with  R  as  follows: 


R2 


oo 


U 

x 

U-  rh 

CO 


(14) 
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Let  us  define" 


f  = 


U 


u 


oo 


Equation  (9)  becomes 


h  = 


1  + 


4f  2\ 


(15) 


(16) 


For  small  values  of  R.r,  the  behavior  of  h  is  parabolic 


E»i--p  (Rjp  (17) 

k  yL 

It  vrill  be  recalled  that  in  Ref.  2  a  Gaussian  distribution  of  the  form 


h  =  e 


(18) 


was  used  as  an  analytic  expression  for  the  initial  enthalpy  profile.  The 
solution  of  the  conduction-controlled  part  of  the  trail  through  an  integral 
method  depends  critically  on  the  accuracy  of  the  profile  used  in  the  neigh¬ 
borhood  of  R  «  0.  Comparing  Eqs.  (17)  and  (l8)  we  find 


o  i  f 


(19) 


o  R=0 


The  factor  f  defined  in  Eq.  (15)  can  be  easily  computed  from  the  fact 
that  the  stagnation  enthalpy  remains  constant  everywhere  in  the  inviscid 


♦Since  f  Ri  1,  the  Howarth  variable  represents  physically  the  shape  of 
the  shock.  Since  f  varies  between  about  0.80  and  1.00,  the  maximum  error  in 
the  approximation  in  Eq.  (l4)  is  about  10  per  cent. 
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flow  field.  Its  value  at  R  =  0  is  given  by 


For  altitudes  between  60,000  and  100,000  ft  and  ve.locities  ranging  between 
15,000  and  25,000  ft/sec,  the  value  of  f  is  almost  a  constant  and  equals  0.77* 
Equation  (20)  is  in  good  agreement  with  the  numerical  findings  of  Ref.  1, 
where  real-gas  tables  have  been  used. 

The  Gaussian  depth  as  defined  in  Eq.  (19)  can  now  be  calculated  with 
k  =  1.27  7l  =  1.18,  and  f  =  0.77-  We  find 

Pxo  =  S  <21> 

Feldman^^  has  chosen  6  ~  1.0  for  all  his  computations,  a  value  found  for 
best  fit  after  fairing  the  enthalpy  curves  from  the  method  of  characteristics 
into  the  Interval  r  <  r^.  His  value  is  well  justified  by  the  above  result, 
since  for  a  hemisphere-cylinder  geometry  in  hypersonic  flow  ~  1.  In 
Fig.  4  several  enthalpy  profiles  have  been  collected,  corresponding  to 
different  flight  conditions  from  Ref.  1.  It  is  seen  that  the  universal 
curve  fits  all  the  numerical  data  well,  especially  for  low  R's  exhibiting 
at  the  same  time  the  correct  trend.*  Feldman's  Gaussian  approximation  is 

*This  curve  is  traced  for  a  value  of  f  at  R  =  0  of  0.77.  The  value  of 
f  at  other  points  was  computed  by  using  the  fact  that  the  stagnation 
enthalpy  remains  constant.  Note  that  for  the  conditions  of  curve  2  in  Fig. 

3  a  value  k  =  1.4  is  appropriate, (1)  corresponding  to  a  Gaussian  depth  of 
1.50,  incompatible  with  curve  2. 
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also  shown,  and  it  is  clear  that  this  approximation  goes  to  zero  faster 
away  from  the  axis. 

The  velocity-deficiency  law  at  the  section  x  =  x^  can  easily  be 
obtained  by  allowing  ourselves  the  assumption  of  constant  stagnation 
enthalpy  in  the  radial  direction.  The  statement  is  an  exact  one  for 
adiabatic  conditions  with  Prandtl  number  equal  to  one. 

We  find 


For  h/hQo»  1  and  U/U^  equaling  a  small  number  compared  to  unity,  the 
above  is  reduced  to 


In  other  words  the  velocity-deficiency  law  is  exactly  the  same  as  the 

enthalpy  law  when  both  are  referred  to  their  value  at  the  cross  section 

x  =  x  and  R  =  0. 
o 

The  method  just  developed  can  easily  be  applied  in  the  case  of  frozen 

flow.  As  was  pointed  out  previously,  at  the  higher  altitudes  the  transition 

point  from  equilibrium  to  frozen  conditions  occurs  closer  to  the  stagnation 

point.  We  may  assume  under  these  conditions  that  the  expansion  of  the  gas 

from  there  on  takes  place  with  a  constant  value  y  calculated  from  either 

r 

Ref.  9  or  Ref.  10.  Assuming  that  in  the  neighborhood  of  the  stagnation 
region  the  shape  of  the  shock  is  not  altered  (the  value  of  k  in  Eq.  (9) 
remains  the  same),  the  Gaussian  depth  p(0)  will  be  given  by  the  relation 
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p(0)  =  J-L-  (24) 

R=0 

(7 ) 

Frozen  conditions  have  a  rather  small  effect  on  the  velocity;  '  therefore 
p(0)  will  increase,  since  7^,  will  be  always  greater  than  7^.  On  the 
other  hand,  the  temperature  at  the  station  x  =  xq  wild  be  lower,  but  because 
of  the  increased  Gaussian  depth  it  will  diminish  in  the  radial  direction  at 
a  slower  rate.  Thus,  all  the  observables  that  depend  on  temperature  will 
probably  remain  high  even  at  large  distances  from  the  axis  of  the  trail. 

The  significance  of  this  behavior  must  not  be  overlooked.  We  have 


fV<°> 

eEtoj 


«  0.85  7p 


(25) 


As  an  example,  consider  the  case  of  flight  at  an  altitude  of  100,000  ft 

and  flight  velocity  of  20,000  ft/sec.  Under  these  conditions  Ref.  11  gives 

a  compressibility  factor  behind  the  normal  shock  of  1.58  and  hence  an  atom 

concentration  of  0-38.  From  Ref.  9  ve  calculate  for  this  concentration  a 

value  of  7  =  1.52,  and  the  same  value  is  obtained  by  using  the  results  of 

r 

Ref.  10.  The  ratio  of  the  two  Gaussian  depths  for  equilibrium  and  frozen 
conditions  is 

PF(0)  _  1.52  _  x  go 

p^oy  -  m b 

A  change  of  approximately  30  per  cent  in  the  degree  of  persistence  in  the 
radial  direction  of  the  high  temperatures  prevailing  at  the  axis  of  the 
trail  could  be  of  importance  for  some  of  the  observables. 

Unpublished  numerical  calculations  for  frozen  flow  by  R.  Goulard 


verify  Eq.  (25). 
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THE  LENGTH  OF  THE  EXPANSION  PART  OF  TIE  TRAIL 

For  the  estimation  of  the  distance  x  one  can  use  the  second-order 

o 

blast-wave  theory  due  to  Sakurai^22^  for  the  pressure  decay  along  the 
distance  x.  We  have 


=  0.1330  +  0.405 

p  x/r 

*  00  '  o 


Equating  the  above  to  one  and  solving  for  x/r  ve  find 


x 

o 

r 

o 


M  2 
00 


The  following  table  compares  the  above  approximation  with  numerical,  results 
obtained  by  the  method  of  characteristics 


Altitude  (ft) 

Velocity  (ft/sec) 

—  =  M  2/ 
r  co/ 4.  5 

Method 

Character: 

AVCO 

60, COO 

17,500 

72 

56 

100,000 

15,000 

53 

4o 

100,000 

20,000 

95 

80 

100,000 

25,000 

145 

100 
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THE  THERMAL  CONDUCTION  PART  OF  THE  TRAIL 

We  proceed  now  to  the  calculation  of  the  conduction-controlled  part 

of  the  trail.  From  now  on  the  measurement  of  distance  x  will  be  made  from 

station  x  =  x  ,  which  in  the  new  notation  will  be  denoted  by  x  =  o. 

0 

Let  us  define  the  length  of  the  trail  L  as  the  distance  from  the  cross 
section  where  the  pressure  has  reached  its  ambient  value  to  a  point  on  the 
axis  where  the  temperature  corresponds  to  a  given  minimum  ionization  level. 
Let  b  be  the  characteristic  radius  of  the  trail  in  the  direction  normal 
to  the  flow.  Assuming  that  the  amount  of  heat  convected  in  the  direction 
x  is  equal  to  the  amount  of  heat  conducted  in  direction  b,  we  may  write 


(AT). 


pUC 


(at). 


(26) 


where  (AT)^  is  the  temperature  difference  over  the  length  L,  (AT)&  is  the 
temperature  difference  between  the  axis  and  the  edge  of  the  trail,  and  U 
is  the  characteristic  velocity  at  the  axis  of  the  trail.  On  the  other  hand, 
the  order  of  magnitude  of  the  drag 


or 


Drag  ■ 


pb2!]2 


2  2 
'Crp  U 
D  o  00  00 


(27) 


where  r)  is  the  body  radius.  After  introduction  of  the  drag  coefficient  in 
the  above  equation,  elimination  of  b  between  Eqs.  (26)  and  (27),  and 
introduction  of  the  Prandtl  number,  we  find 


Or2 
D  o 


00 


u  (AT)t 

00  L 

U 


Pr 

M 


(28) 
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The  characteristic  velocity  U  prevailing  in  the 
con  be  calculated  from  Eq.  (20).  Introduction  of  the 


axis  of  the  trail 
factor  f  yields 


L  Pr 

-  - 

C  r2  M 

Vo 


P  U 

00  OO 


(AT)X 

(atT 


(29) 


The  temperature  ratio  AT^/AT^  depends  solely  on  the  temperature  at  the 
end  of  the  adiabatic  expansion  (a  sole  function  of  the  flight  Mach  number) 
and  the  temperature  corresponding  to  the  lowest  degree  of  ionization  upon 
which  the  length  L  is  based-  Therefore  we  conclude  from  Eq.  (28)  that  for 
a  constant  flight  Mach  number  the  length  of  the  trail  will  be  directly 
proportional  to  the  drag  coefficient,  to  the  cross  section  perpendicular 
to  the  flow,  and  to  the  ambient  density,  which  varies  exponentially  with 
altitude.  It  is  obvious  that  any  detailed  analysis  of  the  problem  can 
yield  no  more  than  a  numerical  constant  in  front  of  the  right  side  of 
Eq.  (28)  and  the  influence  of  the  viscosity  u,  which  is  a  function  of 
temperature.  The  order  of  magnitude  of  the  lengtn  of  the  trail  can  be 
correctly  estimated  from  Eq.  (28). 

We  now  proceed  to  a  more  detailed  analysis  of  the  problem.  The  energy 
equation,  allowing  for  variable  properties  but  neglecting  the  viscous 
dissipation  and  the  components  associated  with  the  radial  velocity,  is 


pU 


dh 

OX 


1  _d_ 
r  dr 


( 


kr 


8T  \ 
or  / 


(50) 
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where  k  is  the  coefficient  of  thermal  conductivity.  Following  Ref.  1  we 


define 


2L  (  — 

r  l  p  U  r  j 
O  \  00  00  o  / 


(31) 


Making  rQ  =  1,  Eq. 


(30)  becomes 


_£ _ !L 

^00  ^00 


Hh/h^ ) 

dX 


1  d 
r  5r 


_1_ 

Pr 


00 


sfh/h^i 


(32) 


In  an  integral  method  we  need  to  satisfy  the  conditions  that  the  total 
energy  crossing  any  station  x  at  all  radii  remains  constant.  Mathematically, 
we  must  have 

00 

Jc  /  r~  “  0 

J  co 
o 

At  very  large  distances  away  from  the  body,  the  velocity  and  enthalpy 
distributions  will  be  Gaussian.  On  the  other  hand,  for  small  distances  the 
profiles  will  be  of  the  form  of  Eq.  (9).  We  choose  the  following  enthalpy 
profile  for  the  study  of  the  conduction-controlled  trail  at  small  x 


a(x) 


OO 


1  ^ 
1  +  T 


3  (x) 


i7^oy 


(34) 


Comparison  of  Eq.  (34)  with  Eq.  (17)  at  x  =  0  and  small  values  of  gives 

4 


7Lk 


P(0)P’(0)  =  ^ 


p(0) 


R-0 


Hence,  in  order  to  approximate  the  initial  enthalpy  profile  we  must  set 
’(0)  =  1. 
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At  x  =  0,  Eq.  (34)  becomes 


,Vi‘ 


1  +  r: 


Equation  (35),  when  plotted  with  f  =  0.77,  7T  =  1.18,  and  k  =  1.27, 

falls  slightly  to  the  right  of  curve  5  and  to  the  left  of  curve  1  in  Fig.  4. 
(it  is  not  shown  there  in  order  to  avoid  crowding  the  figure.)  It  con¬ 
stitutes  a  good  approximation  to  the  initial  profile  throughout  the  whole 
range,  a  property  that  the  Gaussian  distribution  at  the  station  x  =  x 

o 

does  not  seem  to  have.  Equation  (34)  can  now  be  substituted  in  Eq.  (33); 


we  find 


a(x)P*(x)  =  const. 


We  now  take  Eq.  (32)  to  the  limit  for  RR  — 0.  Making  use  of  Eq.  (34),  we 


da(x)  _  _ 4_ 


PrfR=Q  0’(x)P(O)  ^ 


Using  Eq.  (36)  to  eliminate  0*(x),  we  have 


4  a  (x)  u(x,0' 

PrfR=0  a(o)p(o) 


Substituting  Eq.  (19)  into  the  above,  we  find 


da(x) 

dx 


16 -  a2(x) 

PrkSLa(0)  U  a'!  ^ 


The  following  result  emerges  from  Eq.  (38),  recalling  that  ~  k  ~  CR  (Eq.  8). 
The  gradient  of  the  nondimensional  enthalpy,  calculated  at  the  axis 


of  the  trail  at  the  station  x  =  x  ,  depends  only  on  the  drag  coefficient, 
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transport  properties  not  considered.  It  is  worth  noting  that  it  does  not 
depend  on  the  velocity  deficiency  at  R  =  0.* 

Equation  (38)  can  he  integrated  in  a  closed  form  provided  we  make  the 
assumption  of  an  average  constant  value  for  p(x,0)  as  has  already  "been 
pointed  out  in  Ref.  1.  The  result  is 


a(x)  _ 
a(0) 


(39) 


From  Eq.  (37)  it  is  seen  that  the  slope  of  a(x)  is  inversely  proportional 
to  the  Gaussian  depth  fj(0).  For  large  values  of  x,  the  form  of  Eq.  (3*0 
should  become  less  accurate;  a  Gaussian  distribution  will  be  more  appropriate 
there . 

Under  the  assumption  of  constant  thermal  conductivity,  the  heat- 
conduction  equation  (Eq-  (32))  admits  as  an  exact  solution  the  Fourier- 
Bessel  integral 


M  X  ,2 
- —  —  v  x  r  oo 


?  j0(v  "  ^  \  7 

=  i  -  •  x  "  e  ^  J  h(RH,0)  (XRh)  Jo(XRH)d(XRH)  d.\ 

o  o 

Assuming  an  an  initial  condition  h(R^,0)  a  Gaussian  distribution  of 


the  form 


h(RR,0)  -  h(0,0)  e  -RH  /Po 


*The  Gaussian  depth  p(0)  from  Eq.  (19)  was  found  to  be  inversely 
proportional  to  f^_Q. 


Eq.  {ho)  yields  after  integration 


i(O.X) 

i(0,0) 


H 


its: 


^00  fPr^c 


(te) 


^  t'Prp 
00  Kc 


This  is  the  solution  found  by  Feldman  in  Refs.  1  and  2k  by  using  an 
integral  method;  however,  the  above  arguments  show  that  under  the  condition 
Eq.  (4l)  it  is  the  exact  solution. 

If  we  use  Eq.  (35)  in  conjunction  with  Eq.  (21)  as  the  initial  condition 
in  Eq.  (40),  the  first  integral  yields^'^0^ 


00 


1  nl  ... 

j  — ^  trr  wrh^  d 

0  (i  +  4)  D 


(1  + 1  ' 


+  C  ) 

h(0,0)  \  D  K  (\) 

(4  - l) 


(±-l) 

D  '•  r  A- 

*  n 


(43)* 


The  final  result  is 

-hjOj.x) _ L__ 

4(0,0)  /  _1_  ,\ 

>lc"  4  to 


00 
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j  Jo4wVK/l 
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00 


p  Prf^3 
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Unfortunately,  this  last  integral  is  not  readily  available  in  the 
literature  in  a  closed  form,  even  when  0^=1,  and  lengthy  calculations  based 
on  series  expansion  seem  to  be  necessary.  On  the  other  hand,  we  note  that 


*J  and  K,  ,  are  the  Eessel  Functions  of  the  first  kind  of  order 
°  11  -1N1 

IS  }  ,1  , 

zero  and  modified  function  of  the  second  kind  of  order ^  -  1  j  ,  respectively, 


m 
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the  integral  method  based  on  Sq.  (35)  &6  the  boundary  condition  yields  the 

stone  result  at  the  axis  of  the  trail  as  the  exact  solution  based  on  an 

initial  Gaussian  distribution.  This  of  course  results  from  the  parabolic 

character  of  both  distributions  in  the  neighborhood  of  the  axis. 

We  now  turn  our  attention  to  the  possibility  of  obtaining  a  closed-form 

solution  for  the  case  of  variable  thermal  conductivity.  Within  the  range 

of  temperatures  and  pressures  of  interest  here,  the  Prandtl  number  remains 

constant,  whereas  the  viscosity  still  obeys  Sutherland* s  law.*  From 

Fig.  5  it  is  seen  that  the  two  suggested  correlations  of  viscosity  with 

enthalpy  are  in  good  agreement  with  the  exact  values  for  pressures  between 
-1  -2 

10  and  10  atm  for  nondimensional  enthalpies  less  than  100,  whereas  for 

-3 

higher  enthalpies  the  correlations  are  good  for  pressures  between  10  and 
10-4  atm. 

Assuming  in  general  a  power  law  for  viscosity  of  the  form 


p  =  Po  (h/RTo) 


CO 


(45) 


with  Pr  =  0.70,  k  =2.62  and  7^  =  1.20,  we  find 


h(0,x 


h(0,0 


j  -l  +  1-25  “o  (1  +a>)Lh(0,0)j 
^  Poo  U00  [RTo]  ^ 


o  J 


1  +  cu 


(46) 


For  the  calculations  of  Figs.  6-10  and  12,  both  0.25  and  0.28  exponents 
were  used  for  cu  ;  it  was  found  that  the  difference  in  the  results  was 


*For  altitudes  up  to  250,000  ft  and  velocities  up  to  25,000  ft/sec, 
the  temperature  at  the  axis  after  expansion  to  ambient  pressure  is  less 

than  4500°K. 


OIOO'O 


Viscosity  lews  as  power  functions  of  enthalpy 
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small.  In  Fig.  6,  Eq.  (46)  is  plotted  for  a  value  of  to  -0.25 


In  Eq.  (47)  it  is  necessary  to  express  in  grams  per  cubic  cm. 

Figure  7  is  the  result  of  calculations  based  on  Eq.  (47)  for  different 
altitudes  end  velocities  in  terms  of  the  electron  concentration  rather 
than  the  enthalpy. 

We  observe  that  for  the  same  velocity  the  electron  concentrations  are 
higher  for  the  lower  altitudes,  and  at  the  same  time  the  rate  of  decrease 
in  also  slower.  On  the  other  hand,  for  the  same  altitude,  lower  velocities 
correspond  to  lower  electron  concentrations  with  a  faster  rate  of  decay  in 
the  axial  direction.  For  a  given  lower  electron  concentration  the  results 
of  Fig.  7  may  be  used  to  plot  the  length  of  the  trail  versus  altitude  for 
different  velocities.  These  calculations  are  shown  in  Figs.  8  -  10. 


THE  INFLUENCE  OF  THE  TRAILING  SHOCK 

As  was  stated  in  the  introduction,  the  influence  of  the  trailing 
shock  in  the  observables  associated  with  the  conduction-controlled  part  of 
the  trail  will  be  rather  limited.  In  order  to  give  an  estimate  of  its 
importance  in  a  closed  form  we  make  the  extreme  assumption  that  the  state 
of  the  gas  Just  before  the  trailing  shock  is  such  that  the  free-stream 
pressure  has  been  reached.  Some  simple  calculations  entirely  similar  to 
the  ones  that  led  to  Eq.  (9)  yield  the  following  result 


—  Electron-concentration  decay  at  the  axis  of  the  ion 
trail  for  different  altitudes  and  velocities 


concentration  of  10  e/cm 


concentration  of  !09  e/cm 


RM- 2682-1 
50 


From  the  above  we  observe  that  the  equivalent  new  Gaussian  depth  will  be 
smaller  than  the  one  found  when  the  trailing,  shock  was  neglected,  by  the 
factor  (2  -  l/y^)  1  —  0.86.  This  means  that,  apart  from  the  fact  that  the 
temperature  in  the  center  line  at  the  end  of  the  second  expansion  will  be 
higher,  these  higher  temperatures  will  persist  over  a  smaller  range  in  the 
neighborhood  of  the  axis.  Furthermore,  the  temperature  at  the  axis  of  the 
trail  will  decay  faster,  since  we  have  shown  before  that  the  rate  of  decay  ia 
inversely  proportional  to  the  Gaussian  depth. 

Figure  11  shows  the  extent  of  the  influence  of  the  second  shock  under 
the  above  assumptions  for  an  altitude  of  100,000  ft  and  three  different 
velocities.  As  is  evident  on  physical  grounds,  the  second  shock  is 

more  influential  at  the  higher  velocities  where  the  length  of  the  trail 

# 

will  be  somewhat  greater.  In  terms  of  the  nondimen sional  distance  X  , 
as  defined  in  Fig.  11,  for  the  velocity  of  25,000  ft/sec  the  difference  is 
of  the  order  of  17  per  cent,  whereas  for  the  remaining  tvo  velocities  the 
difference  is  negligible.  Recalling  that  the  above  calculations  correspond 
to  rather  extreme  conditions,  we  conclude  that  the  role  of  the  trailing 
shock  in  the  conduction  part  of  the  trail  iB  limited. 

DISCUSSION  OF  RESULTS 

The  present  analysis  has  Bhovn  that  the  determination  of  the  length  of 
the  ion  trail  under  the  assumption  of  thermodynamic  equilibrium  is  controlled 
entirely  by  the  Bhape  of  the  bow  shock  wave  or  the  drag  coefficient  of  the 
object.  For  bluff  bodies  (Cp  of  about  0.$  or  higher)  boundary-layer  effects 
will  be  negligible.  On  the  other  hand,  for  slender,  long  bodies  these 
effects  cannot  be  neglected,  and  therefore  the  present  results  should  be 
understood  as  providing  an  upper  limit  for  the  length  of  the  trail.  It 
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Ghould.  also  be  added  that  for  slender  but  short  bodies  the  influence  of 
the  second  shock  could  be  rather  prominent. 

The  equilibrium  solution  presented  in  this  memorandum  is  certainly 
valid  for  altitudes  less  than  150,000  ft.  For  high  altitudes,  where  the 
assumption  of  frozen  flow  could  be  Justified,  the  basic  form  of  the 
equilibrium  solution  for  the  distribution  of  temperature  or  enthalpy  at 
the  cross  section  at  the  end  of  the  adiabatic  expansion  can  be  maintained, 
based  nevertheless  on  the  frozen  ratio  of  the  specific  heats.  Solution, 
of  the  problem  by  taking  into  account  the  exact  chemistry  of  all  the 
species  present,  including  ambipolar  diffusion  and  electron  attachment  in 
both  the  expansion  and  conduction-controlled  regions,  has  yet  to  be 
obtained. 

From  Figs.  8  -  10  we  observe  that  for  const-ant  velocities  the  trail 
length  is  an  exponential  function  Of  altitude.  As  a  matter  of  fact,  the 
exponent  is  very  nearly  the  same  as  the  one  corresponding  to  the  density 
of  an  exponential  atmosphere.  If  the  altitude  is  denoted  in  thousands  of 
feet,  the  exponent  for  the  density  variation  in  a  simple  isothermal  atmos¬ 
phere  is  about  (-  alt/23.5),  whereas  the  lines  in  Figs.  8-10  correspond 
to  exponentials  equal  to  about  (-  alt/19).  The  steeper  slope  corresponding 
to  the  trail  length  comes  about  from  the  fact  that  for  constant  flight 
velocity  the  value  of  the  enthalpy  needed  for  the  lowest  electron  concen¬ 
tration  upon  which  the  trail  length  is  based  is  higher  for  high  altitudes 
where  the  mass  density  is  low.  Inspection  of  Eq.  (U7)  shows  that  the 
ratio  h(0,0)/h(0,L),  apart  from  its  dependence  on  the  Mach  number  through 
h(0,0),  depends  also  on  the  altitude  through  h(0,L),  as  was  Just  explained^ 
in  this  fashion  the  effective  density  exponent  changes  somewhat. 
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In  Fig.  12  a  plot  is  made  of  the  length  of  the  equilibrium  ion  trail 

for  an  actual  re-entry  monitored  by  three  different  frequencies  corresponding 

7  0  9  3 

to  electron  concentrations  of  10',  10  ,  and  10  electrons/cm  .  As  an 

* 

example,  for  the  last  electron  concentration,  between  the  altitudes  of 

O  _ 1 

100,000  and  60,000  ft,  we  see  from  Fig.  12  a  value  for  (L/C^r^)  —  3000  cm 

2  2 

For  an  object  of  cross  section  r^  —  100  cm  with  a  drag  coefficient  of  one, 
the  length  of  the  trail  will  be  3  km.  For  the  same  object  at  an  altitude 
of  150,000  ft,  the  length  will  be  only  200  m. 

Finally  we  close  with  the  remark  that  for  given  altitude,  velocity, 

2 

observed  length  of  the  ion  trail,  and  the  ballistic  coefficient  (w/C Dr0), 

2 

use  of  a  graph  similar  to  the  one  of  Fig.  12  yields  the  quantity  (cDr0)> 
from  which  the  weight  of  the  re-entering  object  can  be  estimated. 

*10^  electrons/cm^  correspond  roughly  to  a  frequency  of  300  Mcps. 
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